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S enior secondary students cover numerical integration techniques in their 
mathematics courses. In particular, students would be familiar with the 


paper derives these techniques by methods which secondary students may not 
be familiar with and an approach that undergraduate students should be 
familiar with. Secondary students will also find interesting the two-point Gauss 
rule, which is an extension of the trapezoidal rule. There are many applica- 
tions of integral calculus and developing a deeper understanding of some of 
the numerical methods will increase understanding of the techniques. The 
methods chosen in this paper have been investigated as secondary students 
will be familiar with their applications. However, secondary text books and 
teachers may not use the techniques covered in this paper, and this alternate 
approach may increase the understanding of the importance and applica- 
tions of the techniques, as well as increase an appreciation of the beauty of 
mathematics in general. The paper also provides a detailed summary of the 
techniques that will be beneficial for undergraduate students. 

Numerical integration enables approximations to be found for 


where the integral for f(x) cannot be written in terms of elementary func- 
tions. A use of the definite integral is to determine the area between a curve 
and the horizontal axis (see Figure 1 ) . 

In this article the midpoint rule, the elementary trapezoidal rule, the two- 
point Gauss rule and Simpson’s elementary rule will be developed. 


midpoint rule, the elementary trapezoidal rule and Simpson ’s ride. The following 


b 


J f(x)dx 


a 



Z w >/( x *) 
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In this form, A’ is a natural number, w, are called the weights or coefficients and 
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y 



the values of /(x) are called the ordinates. The quadrature rule is not exact for 
every /(x). However, the rule can be exact for some simple functions such as 
1, x, x'. Say we wish to integrate 

i 

J f{x)dx 

o 

by choosing an appropriate xv\ and /(xj) and N= 1, i.e. : 

1 N 

[ /(x)dx = ^a; ] /(x 1 ) 

0 i= 1 

If we consider f(x) = 1, then f(x\) = 1. 

i i 

Now, J f{x)dx = J" 1 dx = 1 

0 0 

So, W\.f{X\) = Wj.l = XV \ 

1 

By equating J f(x)dx and w t .f(x t ), then w 1 = 1. 
o 

Say we now consider /(x) = x, then f(x x ) = x 1 . 


1 1 
J f (x)dx = J x dx = — 

Now, o o 2 

So, ?C|./(X|) = Wj.Xj. 

l ^ j 

By equating J f (x)dx and xv x .j{x x ), then: w x .x x = — . But since xv ] = 1, then x, = — . 


Therefore 


1 

J / (x)c/x = / 




v2y 


which is the midpoint rule on [0, 1]. 


Now, consider the midpoint rule on \a, b\. 

b 

Say, J f {x)dx ~ w x f (xj) 

a 
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If/(x) = 1, then 
Also, 

Therefore, 

Consider /(x) = x, then 


J 1 dx — b — a. 

a 

w l-f(x i) = Wj.l = W 1 

u>i = b - a. 


J x dx = ~° 2 ) 


Now, w 1 .f(x 1 ) = w v Xi, but, Wj = b - a, so, 

[b-a)^ = ^-{b 2 -a 2 ) 

(b — a'j.x j — ^{b — a) (b + a) 


x, = ■ 


b + a 


Therefore, 


D 

J f(x)dx~(b — a)f 


( i \ 

a + b 


v 2 j 


which is the midpoint rule on [ a, b\ . Note that the midpoint rule will integrate 
linear functions exacdy. 


In using the endpoints of the interval when using the simple quadrature rule, 

1 N 

J f{x)dx ~ ^ g ,w i f{x i ) 

0 i= 1 

we can commence with x 1 = 0, x 2 = 1 where N= 2; i.e.: 

i 

J f(x)dx = iv 1 f(0) + w 2 f(V) 
o 

As there are two unknowns, w 1 and w 2 are chosen such that 1 and x can be 
exactly integrated. 

i 

Now, when /(x) = 1 then J f{x)dx-\ 

o 

That is, Wi/(0) + w 2 /(l) = w v \ + w 2 .l = 


as/(0) =/(l)= 1. 
So, Wi = w 2 = 1. 


Letting /(x) = x 


i 

J f(x)dx 


1 

2 


Then, 


w 1 .f( 0) + w 2 .f( 1) = Wi-0 + w 2 . 1 = w 2 


as /(0) = 0 and /(1)= 1. 

Therefore, uk, = — , which leads to w< = — . 

' 2 2 
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Thus 


]f(x)dx = ±[f(0) + /(!)] 


Now, consider die interval \a,b\ where x t = a and x 2 = b. 


Wheny(x) = 1, 

Also, 

Therefore, 

Now, when f(x) = x, 
Also, 

Now, 


u 

J f{x)dx 


— b — a 


iVyf{ a) + w 2 f(b) = Wy + w 2 
Wy + w 2 = b- a. 

b , 

J f(x)dx = —{b 2 -a 2 j 

a 

Wyf{a ) + Wvfib) = w 1 a+ w 2 b. 
Wya + w 2 b = — « 2 ) 

Wi + w 2 = b - a 
Wyd + w 2 b = ^ (b — a^(b + a) 


•••(I) 
• ••(2) 


Multiply Equation (1) by a and subtract from Equation (2): 

b—a . b—a 

Wc, = and w , = 

2 1 2 


Therefore, the general form of the elementary trapezoidal rule on \a,b\ is: 

b j 

\f{ x )dx = ^[f(,a) + f{b)\ 


This Rule will integrate linear functions exactly. 


The two-point Gauss rule is an extension of the trapezoidal rule. N = 2 is 
now considered. This quadrature rule is of the form: 

i 

J f(x)dx = w l f(x l ) + w 2 f(x 2 ) 
o 

and unlike the trapezoidal rule in which x, and x 2 are fixed at the ends of the 
interval, X] and x 2 are not predetermined. 

As there are four unknowns, W\, w 2 , x t and x 2 are chosen such that 1, x, x 2 
and x 3 can be exactly integrated. 

When f(x) = 1: 

i 

J 1 dx = 1 

0 

= W 1 +W 2 
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When /(x) = x. 


1 

J x dx = 
o 


1 

2 


= Wj.Xj +IV 2 .X 2 


When f(x) = x 2 : 


When /(x) = X: 


x 2 dx = — 


= w 1 .x 1 + r</ 2 .x 2 


j" x 3 dx = — 

J 4 


3 3 

= M/j Xj +W 2 .X 2 

In solving these four equations in four unknowns: 

1 


1 1 

x, = 


' ~ 2 2 -ft 


1 1 


x 2 = — + - 


2 2v/3 

This gives the two-point Gauss rule: 

1 ^ (l 1 

J /(x) dx ~ —f 


2 2V3 Z*" 2 2 2x/3 


1 1 

— + - 


The technique to solve the above system of equations is difficult and 
beyond secondary school methods. An appropriate method can be found in 
Kelly (1967, p. 57). 


When considering the general integral: 

b 

J /(x) dx ~ w 1 f(x 1 ) + w 2 f(x 2 ) 

a 

the Two-Point Gauss Rule can be derived similarly by considering the follow- 
ing integrals: 

For /(x) = 1: 

b-a 

Wj + w 2 

b 2 — a 2 

2 

WjXj +w 2 x 2 
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x dx — 


,3 3 

b — a 


— w jXj + w 2 x 2 


b ,4 4 

f 3 b -a 
\ x dx — 

J 4 


3 3 

= WjXj +w 2 x 2 


The solutions to the simultaneous equations are: 

b — a 


b — a 


x, = 


/ - 4 


9 

v 4 

/A ~ \ 

)-a 


-1 

Vi 

1 

vViy 


+ - 


b + a 


+ - 


6 + a 


Hence, 


j f(x)dx ~ 


f 

b-a 

M 

b + a 

b-a , 

( 

b-a 

M 

\ 

b + a 

l 2 

,Vs, 

+ 

2 J 

+ 2 / 

{ 2 

vVs y 

+ 

2 J 


The technique to solve the above systems of equations is difficult, and an 
appropriate method can be found in Kelly (1967, p. 57). 


Simpson’s elementary rule will now be considered, with N= 3. Say on the 
interval [0,1] we choose 

Xj = 0, x 2 = i and x 3 = 1 


for 


1 

J f (x)dx — w,/ (x^ + in,/ (x 2 ) + w 3 f (x 3 ) 


Then 


1 

J /(x) dx = 0 ^/( 0 ) + zc 2 / 


v2y 


+ w 3 /(l) 


Values need to be chosen for w 1 , w 2 and w 3 such that 1, xand x 2 can be exactly 
integrated. 


For /(x) = 1: 


1 

J" f (x)dx = 


= w 1 + w 2 + w 3 


For /(x) = xr. 


i , 

J f{x)dx = — 


2 

1 

= -w 2 +w 3 
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For /(x) = x 2 : 


1 

J f(x)dx 

o 


_ 1 
~ 3 

1 

= — Wa + Wo 

4 2 


Solving the three simultaneous equations: 

1 , 2 
to, = Wo = — and rc 9 = — 
6 '3 


Hence, the elementary Simpson’s rule is given by: 


i j 

J / (x)dx ~ — 


/(0) + 4/ 


V 

v2y 


+ /(D 


Simpson’s rule integrates quadratics as well as cubics exactly. 


If the interval \a,b\ is taken it can be shown that the elementary Simpson’s 
rule is given by: 


0 

J / (x) dx — 


b — a 


f 

( 

\ 

/(«) + 4/ 

a + b 

+m 

V 

K 2 ) 

/ 


A parabola is taken that passes through the following points of the func- 
tion in the interval [ a, b\ : ( a,f( a) ) , ( b,f(b )) and (a+b,f(a+b) ) . The area under 
the parabola then estimates the area under the function. Following is a proof. 


Let F(x) = Ax + Bx+ C be the equation of a parabola, then: 

-2 ■ 


Ax + Bx + C)dx 
Ax 3 Bx 2 


- + - 


- + Cx 


3 2 

A(& 3 -a 3 ) B(b 2 -a 2) J 


+ 


C(b-a ) 


b — a 


b — a 
~6~ 
b-a 
6 

b — a 


[2A(a 2 +ab+b 2 Y?>B(a + b) + bC^ 

(2 Aa 2 + 2 Aab + 2 Ab 2 + ZBa + 3Bb + 6Cj 


Aa~ + Ba + C + Aa 2 +2 Aab + Ab~ +2,Ba + 2 Bb + 4G + Ab~ + Bb + c) 


b — a 


6 


Aa~ + Ba + cj + 4 
/ 

F(a) + 4F " 4 


f ( , \2 / \ A 


A 


a + b 


V 2 j 


+ b 


a + b 


V 2 j 


+C 


+ Ab 2 + Bb+C\ 


f ,\ 
a + b 


v 2 J 


+ F(b) 


To improve the accuracy of applications of the discussed rules, the number 
of sample points can be increased by deriving more complicated rules or by 
dividing the range into many sub-intervals. 
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As can be seen, the numerical integration methods are able to approxi- 
mate a value for a definite integral using the values of the function at points 
within the interval of the integrand. 
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